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We present a current best estimate of the integrated near-infrared (NIR) extragalactic background 
light (EBL) attributable to resolved galaxies in J, H, and K s , Our results for measurements of vlv 

in units of nW irr 2 sr _1 are 1 1 .7+2 6 in 11 - 5 -i 5 in H and 10 -°^o 8 in Ks ~ We derive these new 
limits by combining our deep wide-field NIR photometry from five widely separated fields with 
other studies from the literature to create a galaxy counts sample that is highly complete and has 
good counting statistics out to JHK S ^27 — 28. As part of this effort we present new ultradeep 
K s — band galaxy counts from 22 hours of observations with the Multi Object Infrared Camera and 
Spectrograph (MOIRCS) instrument on the Subaru Telescope. We use this MOIRCS K s — band 
mosaic to estimate the total missing flux from sources beyond our detection limits. Our new lim- 
its to the NIR EBL are in basic agreement with, but 10 — 20% higher than previous estimates, 
bringing them into better agreement with estimates of the total NIR EBL (resolved + unresolved 
sources) obtained from TeV 7— ray opacity measurements and recent direct measurements of the 
total NIR EBL, as well as recent model estimates for the total light from galaxies. We exam- 
ine field to field variations in our photometry to show that the integrated light from galaxies is 
isotropic to within uncertainties, consistent with the expected large-scale isotropy of the EBL. 
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1. Introduction 

The near-infrared (NIR) extragalactic background light (EBL) is the total light from resolved 
and unresolved extragalactic sources in the NIR. This represents the integrated light from all star 
and galaxy formation processes over the history of the universe that has been redshifted into the 
NIR. Some fraction of the NIR EBL can be resolved as the light from individual galaxies (Integrated 
Galaxy Light, IGL), but the fractional contribution from unresolved (and perhaps unresolvable) 
sources is not well constrained. A measurement of the total NIR EBL minus the IGL will provide 
insight into the energy budget of the early universe. 

While direct unresolved measurement of the NIR EBL is technically difficult due to complex 
foregrounds, several authors have reported a measurement of the total (resolved plus unresolved) 
NIR EBL (1; 2; 3; 4; 5; 6; 7; 8) and found it to be a factor of two or more above the IGL obtained 
through source counts (9; 10; 11). This is known as the NIR background excess (NIRBE). The 
spectral energy distribution (SED) of this measured excess, in some cases, appears very similar 
to that of zodiacal light, which suggests there may be a foreground contamination issue. Another 
possible solution to this excess is a large population of undetected faint galaxies and/or population 
III (PopIII) stars contributing a large fraction of the NIR background (see (12; 13) for reviews). 

Since the NIR EBL presents a source of opacity for TeV 7— rays via pair production, the 
density of the background light can, in principle, be measured via direct observation of TeV blazars. 
This method relies on assumptions of the intrinsic blazar spectrum and the SED of the EBL from 
ultraviolet to the NIR, both of which are poorly observationally constrained. Nevertheless, TeV 
7— rays provide an independent estimate of the NIR EBL that can be used to help determine what 
fraction of the background light can be attributed to resolved sources, and how much may arise 
from faint and possibly exotic sources in the early universe. 

2. NIR Background Due to Resolved Galaxies 

The observations and data reduction for the NIR survey data presented in this paper are de- 
scribed in detail in (14), (15) and (16). In Figure [j] we show our completeness-corrected and 
averaged galaxy counts with those drawn from the literature. However, before combining all of 
our data in this way, we first investigated field to field variations for the IGL for our five fields 
over the magnitude range 14.5 < JHK S < 22.5 where all five are highly complete and have good 
counting statistics. We found the IGL over this magnitude range to be consistent across the four 
non-cluster fields (CLASXS, CLANS, CDF-N and SSA13 ~ 7 - 8 nW m~ 2 sr" 1 ) with a 1 a dis- 
persion of ±0.5 nW m~ 2 sr _1 . As such, we find that the IGL is consistent with large-scale istropy, 
an expected signature of the EBL (see (12; 13) for reviews). 

In the A370 cluster field, our IGL measurements were ~ 2 — 3 nW rcT 2 sr _1 higher in all 
bands with the peak contribution to the excess light arising from galaxies at JHK S ~ 17, consistent 
with an excess of L* galaxies at a redshift of ~ 0.4 as in A370. We include the cluster field in our 
average counts for this study because in a survey of a few square degrees such as this, roughly one 
massive cluster will be present. 

Figures [Jd-f) display the same data as in panels (a-c) after dividing by an arbitrarily normal- 
ized Euclidean model (of a = 0.4 in the form N(m) = A x 10 am , and A is a constant) to expand the 
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Figure 1: (a-c) Galaxy counts as a function of apparent magnitude. Average galaxy counts from our 
deep, wide-field data on 5 fields (this study) are denoted by black diamonds. The K s — band counts from 
our Subaru MOIRCS data (this study MOIRCS) are denoted by black squares. Error bars for this work are 
approximately the size of the plot symbols. The counts determined by (14) from the 2MASS (17) field with 
Galactic latitudes of \b\ > 30 are denoted by red asterisks. The counts determined by (14) from the the 
2MASS-6x Lockman Hole survey (18) are denoted by green triangles. The counts determined by (14) from 
three subfields of the UKIDSS (19) LAS are denoted by blue plus symbols. Other data points are taken from 
the studies listed in the plots. The dashed curve shows our error-weighted least squares running average 
(described in Section 2) from which we calculate the NIR IGL. (d-f) The same data as in (a-c) but divided 
by an arbitrarily normalized Euclidean model with slope a = 0.4. 



ordinate and demonstrate where resolved galaxies contribute the most to the IGL. A flat line in (d-f) 
would imply an equal contribution to the IGL at all magnitudes. The areas of positive slope show 
where galaxies contribute a larger fraction to the IGL as one moves toward fainter magnitude. The 
steep negative slope beyond JHK S > 23 — 24 demonstrates the diminishing contribution of resolved 
galaxies to the IGL at the faintest magnitudes. 

The results of our NIR IGL calculation are shown in Figure ^[ Figure ||]b shows the calculated 
slope of the galaxy counts curve as a function of apparent magnitude. Due to simple geometrical 
effects, when the slope drops below a = 0.4 the total light from galaxies begins to converge. 
Figure ||a shows that the vast majority (~ 90%) of the resolved IGL arises from galaxies in the 
range ~ 15 < JHK S < 24. Over this entire range our study is highly complete and has good counting 
statistics. 

Furthermore, from Figure ^]it can be seen that for sources fainter than JHK„ = 28 to make any 
appreciable contribution to the IGL the counts curve would have to rise dramatically over several 
magnitudes beyond the limits of current surveys. Such a rise would quickly result in several (or 
more) galaxies per square arcsecond, rendering any counting exercise impossible due to confusion. 
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Figure 2: (a) Percent contribution to the IGL as a function of apparent magnitude. Solid curves show the 
cumulative contribution from galaxies to the NIR IGL as a function of apparent magnitude. Dashed curves 
show the fractional contribution per each magnitude bin. This demonstrates that the vast majority of the IGL 
(~ 90%) arises from galaxies in the apparent magnitude range ~ 15 < JHK S < 24 where this study is highly 
complete and has good counting statistics, (b) Slope of the galaxy counts curve as a function of apparent 
magnitude. Dash-dotted curves show the measured slope of the galaxy counts curve as a function of apparent 
magnitude. When the slope drops below a — 0.4 the total light from galaxies becomes convergent. 



Faint galaxies certainly exist beyond the limits of the deepest NIR surveys, because at high redshifts 
the faintest apparent magnitudes observed are only probing a few magnitudes fainter than L* down 
the luminosity function. It is unknown whether there is a steep upturn in the luminosity function 
toward faint magnitudes, but if such an upturn exists and faint galaxies contribute significantly to 
the NIR EBL, they would need to be so numerous as to be unresolvable. As such, it may be safe to 
say that the resolvable portion of the NIR EBL (the IGL) has, for the most part, been measured and 
that the most important contribution to the resolved portion comes from galaxies in the magnitude 
range ~ 15 < JHK S < 24, for which the deep wide-field data presented here are optimized. 

3. K s — band Missing Flux 

(20) estimated the missing flux component in the / and //—bands from the faint outer parts of 
galaxies and from galaxies below their detection limits using a histogram of flux in all pixels asso- 
ciated with detected objects. We use a similar method to estimate the flux missed in the K s — band. 
In Figure || we show a histogram of number of pixels versus flux for all pixels associated with de- 
tected galaxies (object pixels) in our MOIRCS K s — band mosaic. Noting the linear trend for fluxes 
~ 0.005 — 0.4 /iJy, we fit a line to the data over this range (blue dashed line). The portion of the 
histogram used in the fit represents 60% of all object pixels in the image and 99% of object pixels 
with fluxes greater than the turnover in the histogram at ~ 0.005 jliJy. The slope of the linear fit 
is -0.86. (20) find a slightly steeper value of -1 by simply estimating the slope by eye. Assuming 
that the true flux distribution for faint pixels (< 0.005 /iJy) continues along the same trend, we 
extrapolate the linear fit to approximate the shape of the histogram when all pixels in the image are 
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Figure 3: Histogram showing Logio of flux in /xJy vs. Logio of the number of pixels at that flux for all 
the pixels associated with galaxies in the Subaru MOIRCS K s — band image. The blue dashed line shows a 
linear fit from ~ 0.005 — 0.4 ji Jy, which includes ~ 60% of all pixels associated with galaxies and ~ 99% 
of such pixels containing fluxes higher than the peak of the histogram. The slope of this line is -0.86. We 
extrapolate the linear fit toward fainter fluxes to the point where all the pixels in the image are accounted for 
(red dashed line). The red hashed area shows the missing flux component corresponding to 1.9 nW m~ 2 sr~' 
in the K s — band. 



accounted for. Using this method, we calculate an estimate for flux missed in the faint outer parts 
of galaxies and in galaxies that are below our detection limits. We find the missing flux component 
to be ~ 22% (1.9 nW m~ 2 sr _1 ) of the total K s — band light from resolved galaxies. 

4. Comparison With Previous Results 

In Figure Q we show our IGL results alongside a summary of measurements of the NIR IGL 
(black symbols) and EBL (red symbols) over the past decade. The purple lines show the modeled 
IGL from galaxies at z < 4 from (21). The blue hashed area of Figure [| shows the allowed (1 a) 
NIR EBL intensity (14 ± 4 nW mT 2 sr _1 ) derived from High Energy Stereoscopic System (HESS) 
observations of TeV blazars (22). (23) use 13 TeV blazars and a grid of NIR background intensities 
to further constrain the NIR EBL and find approximate agreement with the results of (22). As 
noted earlier, however, estimates of the NIR EBL from TeV 7— ray opacity measurements rely on 
assumptions about the intrinsic SEDs of blazars and that of the EBL, both of which are poorly 
constrained observationally. Our results (black diamonds) are some 10 — 20% higher than previous 
estimates of the IGL, which puts them closer to EBL estimates from TeV blazar observations and 
the most recent total NIR EBL measurements. 

A large NIRBE has been found by several groups, with perhaps the most striking result being 
that of (6). This excess, when combined with other NIR and optical background measurements, 
showed an apparent spectral break around 1 /im. The excess was originally attributed to PopIII 
stars, with the break corresponding to the redshifted Lyman limit for these stars at z ~ 10. However, 
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Figure 4: NIR IGL and EBL measurements over the past decade as a function of wavelength (/xm). This 
work (black diamonds) is shown with a summary of measurements of the NIR IGL (other black symbols) 
via integrated galaxy counts (9; 10; 11), model IGL (purple lines) (21), total EBL (red symbols) via total 
light minus stars and zodiacal light (2; 1; 3; 4; 7; 6; 8), and via TeV 7— ray opacity measurements (blue 
hashed area) (22). The range indicated for the 7— ray work shows the 1 a error range. Note that data points 
at 1.25, 1.6, and 2.2 jim have been shifted slightly in their abscissa values for clarity. Arrows are used to 
denote upper and lower limits and otherwise error bars repre sent 1 a confidence levels. Our results bring the 
measurement of the IGL into better agreement with TeV 7— ray observations and the most recent total NIR 
EBL measurements. The lower limits on our data points show the 1 a error estimates associated with the 
galaxy counts integration described in Section 2, while the upper limits show these same 1 a error estimates 
plus the missing flux component derived in Section 3 



a search for the possible absorption imprint of this break on the 7— ray SED of blazars did not find 
evidence for such a feature (24). 

5. Summary 

Our new results for the NIR IGL place the best current constraints on the total NIR light from 
resolved galaxies and serve as a new lower limit to the total NIR EBL. While these results are in rel- 
ative agreement with previous measurements, our numbers are 10 — 20% higher, bringing them into 
better agreement with those derived from 7— ray experiments and the most recent measurements of 
the total NIR EBL. 

We find the IGL to be roughly isotropic, consistent with the expectation of large-scale isotropy 
in the EBL. We confirm that the starlight subtraction for the most recent total NIR EBL measure- 
ments is correct (the reader is referred to Keenan et al. (25) for a detailed description of this 
measurement). 

While our measurements cannot rule out the existence of a NIRBE due to PopIII stars or other 
exotic early universe objects, our new lower limits on the IGL and the upper limits found from TeV 
7— ray experiments (22) could now be considered in rough agreement with the most recent total 
NIR EBL measurements in the /—band, and in near agreement in the K— band. 



6 



The Resolved Near-Infrared Extragalactic Background 



Ryan Keenan 



References 

[1] M. G. Hauser et al., The COBE Diffuse Infrared Background Experiment Search for the Cos- 
mic Infrared Background. I. Limits and Detections, 1998, ApJ, 508, 25 

[2] E. Dwek & R. G. Arendt, A Tentative Detection of the Cosmic Infrared Background at 3.5 \lm 
from COBE/DIRBE Observations, 1998, ApJL, 508, L9 

[3] V. Gorjian, E. L. Wright, & R. R. Chary, Tentative Detection of the Cosmic Infrared Back- 
ground at 2.2 and 3.5 Microns Using Ground-based and Space-based Observations, 2000, 
ApJ, 536, 550 

[4] E. L. Wright, DIRBE minus 2MASS: Confirming the Cosmic Infrared Background at 2.2 Mi- 
crons, 2001, ApJ, 553, 538 

[5] T. Matsumoto, Observations of the Near Infrared Extragalacric Background Light, 2001, in 
IAU Symposium, Vol. 204, 87 

[6] T. Matsumoto et al., Infrared Telescope in Space Observations of the Near-Infrared Extra- 
galactic Background Light, 2005, ApJ, 626, 31 

[7] L. Cambresy et al., The Cosmic Infrared Background at 1.25 and 2.2 Microns Using DIRBE 
and 2MASS: A Contribution Not Due to Galaxies?, 2001, ApJ, 555, 563 

[8] L. R. Levenson, E. L. Wright & B. D. Johnson, DIRBE Minus 2MASS: Confirming the CIRB 
in 40 New Regions at 2.2 and 3.5 \im, 2007, ApJ, 666, 34 

[9] R Madau & L. Pozzetti, Deep galaxy counts, extragalactic background light and the stellar 
baryon budget, 2000, MNRAS, 312, L9 

[10] T. Totani et al., Diffuse Extragalactic Background Light versus Deep Galaxy Counts in the 
Subaru Deep Field: Missing Light in the Universe?, 2001, ApJL, 550, L137 

[1 1] R. I. Thompson, Star Formation History and Other Properties of the Northern Hubble Deep 
Field, 2003, ApJ, 596, 748 

[12] A. Kashlinsky, Cosmic infrared background and early galaxy evolution, 2005, Physical Re- 
ports, 409, 361 

[13] M. G. Hauser & E. Dwek, The Cosmic Infrared Background: Measurements and Implications, 
2001, ARA&A, 39, 249 

[14] R. C. Keenan , L. Trouille, A. J. Barger, L. L. Cowie, & W. H. Wang, An Extremely Deep, 
Wide-Field Near-Infrared Survey: Bright Galaxy Counts and Local Large Scale Structure, 
2010, ApJS, 186, 94 

[15] W. H. Wang, A. J. Barger & L. L. Cowie, Ultradeep Near-Infrared Observations of Goods 
850-5, 2009, ApJ, 690,319 



7 



The Resolved Near-Infrared Extragalactic Background 



Ryan Keenan 



[16] A. J. Barger, L. L. Cowie, & W. H Wang., A Highly Complete Spectroscopic Survey of the 
GOODS-N Field, 2008, ApJ, 689, 687 

[17] M. F. Skrutskie et al., The Two Micron All Sky Survey (2MASS), 2006, AJ, 131, 1163 

[18] C. A. Beichman et al., A Deep 2MASS Survey of the Lockman Hole, 2003, AJ, 125, 2521 

[19] A. Lawrence et al., The UKIRT Infrared Deep Sky Survey (UKIDSS), 2007, MNRAS, 379, 
1599 

[20] R. I. Thompson et al, Evidence for a z < 8 Origin of the Source-subtracted Near-Infrared 
Background, 2007, ApJ, 657, 669 

[21] A. Dommguez, et al., Extragalactic background light inferred from AEGIS galaxy -SED -type 
fractions, 2011, MNRAS, 410, 2556 

[22] F. Aharonian et al., A low level of extragalactic background light as revealed by y-rays from 
blazars, 2006, Nature, 440, 1018 

[23] D. Mazin & M. Raue, New limits on the density of the extragalactic background light in the 
optical to the far infrared from the spectra of all known TeV blazars, 2007, A& A, 471, 439 

[24] E. Dwek, F. Krennrich, & R. G. Arendt, Is There an Imprint of Primordial Stars in the TeV 
y-Ray Spectrum of Blazars?, 2005, ApJ, 634, 155 

[25] R. C. Keenan , A. J. Barger, L. L. Cowie, & W. H. Wang, The Resolved Near-Infrared Extra- 
galactic Background, 2010, ApJ, 723, 40 



8 



